Actinide Topological Insulator Materials with Strong Interaction 
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Topological band insulators have recently been discovered in spin-orbit coupled two and three 
dimensional systems. In this work, we theoretically predict a class of topological Mott insulators 
where interaction effects play a dominant role. In actinide elements, simple rocksalt compounds 
formed by Pu and Am lie on the boundary of metal to insulator transition. We show that interaction 
drives a quantum phase transition to a topological Mott insulator phase with a single Dirac cone 
on the surface. 



Topological insulators are new states of quantum mat- 
ter with a full insulating gap in the bulk and gapless 
Dirac fermion states on the surface[ll-Q. Spin-orbit cou- 
pling plays a dominant role in topological band insulators 
HgTe, Bii_a;Sb:E and Bi2Te3 class of materials. The ba- 
sic mechanism for all topological insulators discovered so 
far is the band inversion driven by spin-orbit couplingQ. 
Soon after the discovery of topological band insulators, 
the concept of a topological Mott insulator has been pro- 
posed theoretically^, where interaction effects are re- 
sponsible for topological insulating behavior. The inter- 
play between electron interaction and spin-orbit coupling 
could enhance the bulk insulating gap, which is impor- 
tant for device applications. Magnetic and superconduct- 
ing order driven by interactions could co-exist with topo- 
logical order, and give rise to exotic excitations such as 
magnetic monopoleQ, Majorana fermion and dynamic 
axion field Q. 

In this work we theoretically predict a new class of 
topological Mott insulators in actinide compounds with 
electrons in the 5f valence orbitals. Interaction and 
spin-orbit coupling energy scales are comparable in these 
compounds, making them ideal material systems to in- 
vestigate the interplay between these two important ef- 
fects. Simple rocksalt compounds formed by Pu and 
Am lie on the boundary of metal to insulator transition, 
where the 5f electrons change from itinerant to local- 
ized characters 0. Using AmN as a example, we show 
that this rocksalt compound is metallic without interac- 
tions, but an insulating gap is induced by interactions. A 
band inversion between the 6d and the 5f drives the sys- 
tem into a topological insulator phase with a single Dirac 
cone on the surface. The predicted actinide topological 
insulators are binary compounds with a simple rocksalt 
crystal structure and reasonably large energy gaps, pro- 
viding an ideal platform to investigate interacting topo- 
logical insulators. Standard topological band insulators 
in the Bi2Te3 class of materials have bulk carriers in- 
duced by disorder. In contrast, the actinide compounds 
predicted in this work have stronger three dimensional 
ionic bonds. As a result, they may have fewer bulk de- 



fects; in fact, earlier photo-emission experiments indicate 
vanishing density of states at the fermi level [Toj. 

Americium and plutonium are actinide elements with 
electron configurations of S/^Gd^Ts^ and 5f^6dP7s^ . 
Many of the Am mono-pnictides AmX (X=N, P, As, 
Sb, and Bi) and Pu mono-chalcogenides PuY (Y=Se 
and Te) are known to be nonmagnetic and narrow- 
gap semiconductors [lol - fisj . At ambient conditions, 
they crystallize in the rocksalt structure. The rock- 
salt structure has face-centered cubic (FCC) symmetry 
with FmSm space group, shown in Fig. lA with AmN 
as an example. It can be represented as two inter- 
penetrating FCC lattices with a two-atom basis. The 
first atom (Am) is located at each lattice point, and 
the second atom (N) is located half way between lat- 
tice points along the FCC unit cell edge. The bravais 
lattice vectors are ai = {0, a/2, a/2), 02 — {a/2,0, a/2) 
and Us = {a/2, a/2,0). Its Brillouin zone (BZ) has three 
independent time-reversal-invariant-momentum (TRIM) 
points, which are r(OOO), X(7r7r0) and L(7r00) as shown 
in Fig. IB. 
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FIG. 1. Crystal structure and brillouin zone. (A) Rocksalt 
crystal structure of AmN with Fm3m space group. Am and 
N separately form a face-centered cubic(fcc) sub-lattice, and 
these two sub-lattices are at a distance of a/2 away with each 
other along one cubic edge. (B) Brillouin zone for AmN of 
rocksalt structure. It has three independent time-reversal- 
invariant points are r(OOO), L{nOO), and X(7r7rO). 
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In order to study the electronic structure of these ma- 
terials, ab-initio calculations are carried out in the frame- 
work of local density functional approximation(LDA) 
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supported by ESTATE package [l^ within the plane- 
wave pseudo-potential method. The ultra-soft pseudo- 
potential[l3l is employed for localized orbitals, i.e., 
Am(Pu) 5f and N 2p orbitals. The k-point grid is taken 
as 16x16x16 and the kinetic energy cutoff is fixed to 
340eV for all the self-consistent calculations. The spin- 
orbit coupling (SOC) interaction is taken into account 
because both Am and Pu have very heavy nuclei. Due 
to the localization of 5f electrons, LDA+U method [Tsj 
is used to deal with strong correlations in these com- 
pounds with U = 2.beV for Am [ll| and J7 = i.QeV for 
Pu[l3|. We use the experimental lattice constants in our 
calculations [191. 
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Because all the AmX and PuY have similar electronic 
structures, here AmN is taken as an example. The band 
structures for AmN with U = QeV and U — 2.5eV are 
shown in Fig. 2A-B separately. The strong interaction U 
is responsible for the metal to insulator transition. If we 
ignore the effect of U in the calculation, there are residual 
5f states at the Fermi level {Ep), indicating metallic be- 
havior (Fig. 2A). Only by taking into count oiU = 2.5eV, 
those residual 5f states are lifted away from Ep, and we 
can correctly predict AmN to be a insulator (Fig. 2 B), 
consistent with the experimental finding [l^. We will 
elaborate this point in the following energy band anal- 
ysis. First of all, we consider the Am-X (Pu-Y) bonding. 
One / and two s electrons of Am bond with the three p 
electrons of pnictides, while two s electrons of Pu bond 
with the four p electrons of chalcogenides. As a result, 
there arc six / electrons left in the ground states for 
both AmX and PuY before we consider the hybridiza- 
tion between 5/ and 6d states in Am or Pu. Since these 
strong bonding states are further away from Ep than 
the low energy Am-Am and Pu-Pu bonding states, we 
ignore them in our band structure analysis in Fig. 3. In 
the following we will just consider the Am-Am and Pu- 
Pu d and / low energy bonding states. The 6d orbital 
first splits into eg and t2g states due to the large crystal- 
field splitting (CFS) and then splits again into F^ and 
F^ subbands due to a small spin-orbit coupling (SOC) 
(Fig. 3)^. The f states split into S/^/^ and 5/^/^ grates 
with angular momentum of 5/2 and 7/2, due to the SOC 
interaction. For the U = 2.5eV case, the Coulomb U 
additionally shifts the 5/^/^ states away from i?_F[ll|, 
leading to the metal to insulator transition (Fig. 2A-B). 
Therefore, only the 5/^^^ states are fully occupied with 
the six / electrons with a non-magnetic ground state, 
before hybridizing with d orbital. CFS further splits the 



FIG. 2. Band structure and topological surface states. Band 
structures for AmN at (A) U = OeV and (B) U = 2.5eV. 
The line width in these fat band structures corresponds to the 
projected weight of Am's d orbitals on Bloch states. The f 
states are divided into two parts with angular momentum 5/2 
and 7/2 due to the spin-orbit-coupling interaction. The part 
with angular momentum 5/2 is occupied, and the other part 
with angular momentum 7/2 is partially occupied with U — 
OeV. For U — 2.5eV case, these partially occupied f states 
of angular momentum 7/2 move further toward high energy 
and become completely unoccupied. (C) Energy spectrum of 
the ab-initio calculation based TB hamiltonian with the open 
boundary condition along '+z' direction. One surface Dirac 
cone is clearly seen at F point as red lines dispersing in the 
bulk gap on the AmN [001] surface. In the inset, the projected 
two dimension BZ is shown. 




FIG. 3. Band sequence. Schematic diagram of the evolution 
from the atomic d;f orbitals of Am or Pu into the conduc- 
tion and valence bands of AmX or PuY at the F point. The 
black dashed arrow represent the effect of turning on chem- 
ical bonding, CFS, SOC and U (see text). The red dashed 
line represents the Fermi energy. The low energy bands in the 
green box explains the topological physics. 

At Ep, the most important feature is the downward 
dispersion of the F^ suband of Am's d band, it crosses 
the F^ subband of the 5/^/^ bands along T — X direc- 
tion, leading to a band inversion near the X point in 
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TABLE I. Energy gap and parity eigenvalue. Eg indi- 
cates the energy gap. F, X and L are the time-reversal- 
invariant (TRIM) points in Brillouin zone. ' — ' and '+' in- 
dicate the 'odd' and 'even' parity product for all the occupied 
bands. 'Tot.' means the total parity product of all TRIM 
points. 



the BZ (Fig. 3). A energy gap is formed at the cross- 
ing point due to the f-d hybridization. Because the d 
and f orbitals have 'even' and 'odd' parities at three X 
points, this inversion of bands with opposite parity signs 
determines the non-trivial topology of these materials. 
This d-f band inversion phenomenon is a generalization 
of the band inversion between the s and p bands in HgTe 
quantum wells Q and between p bands of opposite par- 
ities in Bi2Te3 class of materials [22|. We calculate the 
parity at all TRIM points [H (Table 1) and predict that 
AmN and all other AmX materials in this family are non- 
trivial topological insulators with their band gaps shown 
in Table 1. The LDA+U method cannot obtain a energy 
gap for PuTe, which was measured experimentally to be 
0.2eV at ambient pressure and OAeV with pressure up to 
5 GPa[l3,[l3|- Suzuki and Oppeneer employ dynamical 
mean-field theory (DMFT) method and successfully re- 
produce the energy gap for PuSe and PuTep^.[lH. With 
the same low energy physics as AmN, we can conclude 
that PuSe and PuTe are non-trivial TIs. 

According to the ab-initio calculation, the two low en- 
ergy bands can capture the salient topological fea- 
tures of these materials. Therefore, we propose the fol- 
lowing 8x8 Hamiltonian: 



H 



Hd Hdf 



(1) 



Hd and Hf are 4x4 Hamiltonian for Pj and P^ repre- 
sentations with an effective angular momentum j = 3/2. 
We first write down the continuum Hamiltonian near the 
P point of the BZ and then extend it to a tight binding 
(TB) model for the full BZ. The effective Hamiltonian 
for P^ bands is the Luttinger model (2^ : 



Adif),Bd(f),Cd{f),Dd{f) are constants and J3/2 are spin 
matrices for angular momentum 3/2. The first three 
terms have continuous rotational symmetry; the last term 
breaks continuous rotational symmetry while preserving 
cubic symmetry, which is responsible for the band inver- 
sion only along T — X direction rather than other direc- 
tions in the BZ[2^. The off-diagonal Hamiltonian Hdf 
connects d and / bands of opposite parities, therefore, 
the first order term can be written as Hdf — F{k ■ J3/2) 
with continuous rotational symmetry, F again is a pa- 
rameter. 

Because the 5/ and 6d orbitals are localized, wc con- 
sider only the nearest neighborhood hopping of Am- Am 
and Pu-Pu atoms. Each Am (Pu) atom has twelve near- 
est neighbors located at ±ai, ±a2, ±a3, ±(a2 — ai), 
±(02 — 03) and ±(a3 — ai). So we can extend the above 
continuum model to the TB lattice model and reproduce 
the low energy bands of these materials well, with AmN 
as an example [25t. Besides this single particle part, our 
Hamiltonian also includes the on-site interaction param- 
eter U for electrons on the f orbital. One important ef- 
fect the U interaction is to renormalize the parameters of 
the single particle Hamiltonian, and to induce a metal- 
insulator transition. On the insulator side, we can effec- 
tively fit our model with renormalized parameters to the 
first-principle calculations. 

To observe the surface states, we construct a ab-initio 
calculation for the semi-infinite AmN system with open 
boundary condition in '+z' direction on the basis of 
maximally localized Wannier functions (MLWFs) f26l. [27j . 
We can obtain the surface Green's function of the semi- 
infinite system itcratively. Finally the local density of 
states(LDOS) can be calculated with the imaginary part 
of these Green's functions. From these LDOS, we can ob- 
tain the dispersion of the surface states, shown in Fig. 2C. 
A single Dirac cone at f clearly shows up in the en- 
ergy gap, which demonstrates the non-trivial Z2 topology 
of the material in addition to the parity analysis. The 
Fermi-velocity for the Dirac cone is vp — 1-3 x lO^m/s, 
a little smaller than Bi2Se3[23. To confirm this result, 
we also take open boundary condition in z direction in 
the analytical TB model (1)[1^, from which we have ob- 
served the same single surface state Dirac cone at P point 
of the BZ as in the ab-initio calculation based TB hamil- 
tonian. 



Hd{f) = Cd(f) 



Bdif)ik ■ J3/2f (2) 
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